On a new approach to optical solitons in dielectric fibers 
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in ■ Abstract 

In this letter we describe a new and simple approach for modulating the elec- 

£h ] tromagnetic wave propagating in a dielectric medium in the form of solitons by 

considering the torque developed between the induced dipoles in the medium 

and external field without taking into account the nonlinear Kerr effect. A 

reductive perturbation method deduces the Maxwell equation coupled with 

o 

£ — ' the Larmor equation of the torque to the derivative nonlinear Schrodinger 

o" 

equation that admits optical solitons in the medium. 
PACS number(s): 42.65, 42.65. T, 42.81, 42.79.S 

•i-H ' 

> 

•i-H 

X 



Typeset using REVTeX 



"nld@bdu.ernet.in 



Hasegawa and Tappert |jj showed that under slowly varying amplitude electromagnetic 
(EM)-pulse propagating in a nonlinear Kerr fiber medium is governed by the completely 
integrable nonlinear Schrodinger (NLS) equation %A X +A U +/J*\ A \ 2 A = that admits N-soliton 
solutions (/x=constant). This was derived from Maxwell equations under the assumption of 
weak linear dispersion and A is the slowly varying amplitude of the electric field of the 
electromagnetic wave (EMW). Later this was experimentally verified by Mollenauer et al 
||. For more channel handling capacity, it is necessary to transmit pulses at the order 
of sub-picosecond and femto-second frequency levels. But the propagation of such ultra 
short pulses (USP) experience higher order effects like third order dispersion (TOD), self 
steepening (SS) and stimulated Raman scattering (SRS). In this case the wave propagation 
is described by the higher order nonlinear Schrodinger (HNLS) equation iA x + A tt +fi\A\ 2 A+ 
i(^iA tu +fi 2 \A\ 2 A t + /! 3 |A|^A) = 0. For example, when fi = and when the two inertial 
contributions of the nonlinear polarization namely, the stimulated Raman scattering (SRS) 
and self steepening (SS) are equal (/i 2 = /z 3 = 1) in the absence of third order dispersion 
(TOD) (y^! = 0), the HNLS reduces to the completely integrable soliton possessing derivative 
nonlinear Schrodinger (DNLS) equation || iA x +A tt +i(\A\ 2 A) t = 0, and, however, when/i ^ 
0, it reduces to the completely integrable mixed derivative nonlinear Schrodinger (MDNLS) 
equation [|[ %A X + A u +^\A\ 2 A +i(\A\ 2 A) t = that also admits solitons. When both the 
TOD and SS effects are included one obtains the Hirota equation for specific parametric 
choices (//i = l,/i 2 = ±6,/i 3 = 0) %A X + A tt +fi\A\ 2 A + i(A ttt ±6\A\ 2 A t ) = which 
explains the propagation of ultra short pulses in the form of soliton. Later Sasa and Satsuma 
proved that for a particular choice of parametric relations for TOD, SS and SRS (fii = 
l,//2 = 6,/i3 = 3), the soliton propagation is supported by another completely integrable 
HNLS equation %A X + A tt +{i\A\ 2 A + i(A m +Q\A\ 2 A t + 3\A\ 2 A) = 0. 

The propagation of optical pulses in birefringent fibers have become very useful in 
the context of nonlinear directional couplers and a lot of work has been carried out re- 
cently, in this direction where the dynamical equations governing the propagation of signals 
in the form of optical solitons reduce to the two coupled nonlinear Schrodinger (CNLS) 



family of equations |]|| iA jx + A jtt +/i[ELi l^fc| 2 ]^i +i\piAj XXX + jU 2 (Efc=i I^M 2 )^j* 
+A i 3(5Zl=i( 1^4-A: | 2 )rc^4j] = 0, where j = 1,2. This has been further extended to iV-signals 
H . Very recently, the study of propagation in birefringent optical fibers introduced the new 
concept of shape changing solitons that share energy amongst themselves during propagation 
|0||- This energy switching behaviour of optical solitons has been used for constructing all 
optical logic gates |L(J • The above models (both single and coupled) support propagation of 



optical pulses in the pico- and femto-second ranges that emerge from high intensity lasers. 
Now it has also been experimentally proved that visible white light emerging from an incan- 
descent source, namely a quartz-tungsten-halogen bulb when propagating through the photo 
refractive crystal Sr / r5 Ba 2nNb 2 OQ, optical solitons are formed by self trapping [nj. In an- 



other direction self modulation of quasi-monochromatic EMW into spatially coherent optical 
solitons in a dielectric medium is described by the NLS equation iA t + A xx + 2\A\ 2 A = 
and their higher order and coupled versions, however with an interchange of the time(t) and 
space(x) variables in their derivatives [[L^,[nj. 

In this letter we describe a new and simple approach to modulate the EMW propagating 
in a dielectric medium in the form of spatially coherent optical solitons without taking 
into account the nonlinear Kerr effect in the medium as has been done so far but only by 
considering the torque associated between the induced dipole moment of the medium and 
the electric field component of the EM field in the form of Larmor equation. When the 
dielectric medium is exposed to the weak EM field the induced dipoles, due to polarization 
when continuously subjected to the external field, give rise to torque. To the lowest order 
of perturbation the torque due to the electric field component of the EM field is due to the 
interaction of it with the induced electric dipole moment of the medium. When the medium 
is anisotropically polarizable the induced dipole moment p and the electric field E are not 
parallel in general and the torque is given by T = p A E. The angular momentum vector 
is the only physical quantity that can define a unique direction for the polarised atom of 
the dielectric medium. The classical time evolution equation for the angular momentum 
can be constructed by writing the time derivative of the angular momentum as equal to 



the above torque. Thus we have the gyroscopic equation or the Larmor equation for the 
electric field in the form — -j^- = 7P A E |TJ| , where the polarization P is the sum J2i P« 



over all the atoms in a unit volume and 7 is the Larmor frequency. Even though existence of 
anisotropic polarizability in the dielectric medium is necessary for the presence of an electric 
field induced torque, it is not sufficient and therefore in addition, it is required that the 
angular momentum vector should be neither parallel nor orthogonal to the electric field. 

The propagation of EMW in the dielectric, when there are no free stationary and moving 
electric charges, is described by the Maxwell equations : 

V-D = 0, V-B = 0, (la) 

™ B = -£- v ^f- « lb » 

where the fields E, B, H and D, which are represented by 3-component vectors, bear the 
usual text book meaning and e is the dielectric constant of the medium. As we consider a 
dielectric medium for EMW propagation here the constitutive relations between B and H 
can be written as H = — , where /1 is the magnetic permeability of the medium. Now taking 
curl on both sides of the first of Eq. ([[]]) and using the second of Eq. ([TE]) , and after using the 
relations for H and D = eE + P, we finally obtain J^[E + ±P] = c 2 [V 2 E -V(V • E)], where 
c = -j= and P = (P x ,P y ,P z ). To understand the nature of propagation of EMW in the 
dielectric medium we now have to solve Maxwell equations coupled with Larmor equation. 
As fiber is the convenient form of dielectric medium for propagation, it is more meaningful 
to study the one dimensional version of these equations (say along ^-direction) . We further 
consider x as the direction of propagation of EMW, which is taken parallel to the length of 
the dielectric fiber. Thus we have the following set of coupled equations to be solved. 
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^j[E(i,() + ^P(i,()]=c 2 



dx 2 dx 2 



(2a) 



and 

dP(x,t) 



P(i,t)AE(i,t), (2b) 



dt 

where e = (1, 0, 0) and time t has been rescaled and velocity c redefined while writing Eq.( 



Now, we solve the set of coupled equations (Q) using a reductive perturbation method 
originally developed by Taniuti and Yajima [|15j to understand the modulation of the slowly 
varying EMW due to nonlinear interactions in the medium. When the amplitude varies 
slowly over the period of the oscillation, a stretching transformation allows us to separate 
the system into a rapidly varying part associated with the oscillation and a slowly varying 
one such as the amplitude. Nonlinear terms modulate the slowly varying amplitude of the 
wave. Then a formal solution is given in an asymptotic expansion about a uniform value. 
A similar study has been performed in the context of EMW propagation in ferromagnetic 
medium and different types of EM solitons were obtained depending on the nature of the 
medium [0,0. We assume that for the set of coupled Eqs.(H) there exists a formal solution 



expanded asymptotically in terms of a small parameter e in the neighbourhood of a constant 
value. For this we expand the components of the electric field and polarization of the medium 
in terms of a small parameter e as ^ x (^,t) = JF +e!Ff +e 2 J r $+ ..., jF Q (^,r) = e^ [jFf 
+£^2 +•••]) where a = y,z. Here T stands for the electric field E(x,t) as well as for the 
polarization vector P(x, t). Due to the anisotropic character of the molecules in the medium 
as originally assumed here we have expanded the electric field and the polarization vectors in 
a nonuniform way about the constant values Eq and Pq respectively. The expanded electric 
field and polarization are functions of the slow variables r and £ introduced through the 
stretching of time (r = e 2 t) and the space variable (C, = e(x — vt)), to take care of the slow 
variation in amplitude ||18|| . Here e is the same perturbation parameter introduced earlier 
and v is the group velocity of the propagating EMW. 

We now substitute the expansions of E and P in the component forms of Eqs.(|2|) and 
collect terms proportional to different powers of e and try to solve the resultant equations. 
On solving the resultant equations at 0(e°), we obtain E = — ^ and Ef = kP", where 
k = (E /P ) = e , c 2_ v 2y Then solving them at 0{e l ), after using the results of 0(e°) we 



finally obtain E'f = — - and 
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2v(2c 2 -v 2 ) 
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Now, we define 



^={E\-iE{), 



i> 



El 



(4) 



The above definitions suggest that Ef — E\ + E\ , which is in accordance with the nonuni- 

after a single differentiation, can be combined together 



form expansion. Eqs.(3b) and 

to give after some algebra and rescaling the variable r the following single equation in terms 

of the newly defined variable ip, 



iip T + ip^ + i\[\ ip \ 2 ip]t 



0. 



(5) 



where A = — +e - — -. It may be verified that on using the definitions (|J), Eq.(|3*a]) can 
be written in the form of Eq.(|5]). Eq.(0) is the completely integrable derivative nonlinear 
Schrodinger (DNLS) equation which admits N-soliton solutions. The multi-soliton solutions 
of Eq. (|5|) can be obtained in the framework of the inverse scattering transform (1ST) method 
M. It is moreover easy and straightforward to construct the N-solitons using Hirota's bi- 
linearisation procedure JTTJ. Following ref. [EO], we consider the nonlinear transformation 



v> 



9f* 



so that Eq. (J5|) can be written in terms of the following bilinear equations. 

(tD T + D\)g.f = 0, 
(g.f)(iD r + Dl)+2D^g.f)D^\ (/*./) 
-2(g.f*)DJ(f.f) + i\ [3Ds(g.f) + gfD^} (g*.g) = 0, 



(6a) 



(6b) 



where D? D% (g ■ f) = |£ - £| [| - ^j g(£, r) ■ /(£', r')| a t e=e , T>=T are Hirota's 
bilinear operators |19|| . In order to obtain soliton solutions, we assume the series expansions 

for^,/andras^ = x^i + X 3 ^3 + -, / = l + X 2 / 2 + X 4 /4 + ... and/* = l + x 2 /2+X 4 / 4 * + - 
where \ is an arbitrary parameter. Now, for example to construct a one-soliton solution we 



set g = X9i, / = 1 + X 2 f2 and /* = 1 + x 2 fl> an d then collect terms with similar powers of 
X and solve the resultant equations which admit the following solutions for g, f and /*. 

g(£,r) = exp[rji]. (7a) 

f ^- 1 + W^W exp[m + " l] ' (7b) 

where 771 = 771^ + ir\\i = K\T + f2i£ + 77} , K 1 = iVL\ and Vti = &x R + i£lu, and the constant 

3l 

f* 

ip = Qsech(rj 1R + A )tanh(r] 1R + A ) , (8) 



^1 = Vir + ivii ■ Thus the one soliton solution for ip = ^ is found to be 



where Q = -^-exp[ir)u + A], A = r][J, + A and A = |^(|^-)- Similarly we can find two, 
three ... and N-soliton solutions. As the details of constructing multi-soliton solutions are 
very lengthy and the form of iV-soliton solution is very cumbersome we are not presenting 



the details here and interested readers can refer to ref. |2(J. Using the one soliton solution 
given in Eq.(|) in the definition of ip (Eqs.(fD) the components of the electric field vector can 
be constructed. Thus the optical one soliton in terms of the electric field vector component 
of the EM field can be written as 

— P x 1 

El = = -exp(2A)sech 2 (r] 1R + A )tanh 2 \r] lR + A ), (9a) 

E? = kP? = --exp(A)[Re(Q)5 xa - Um(Q)5 ya ]sech(r] 1R + A )tanh(r] 1R + A ), (9b) 

where Re(Q) and Im(Q) represent the real and imaginary parts of Q respectively and S xa 
and S ya are Kronecker delta functions. 

Eqs.(^) show that the electric field component of the propagating plane EMW is mod- 
ulated in the form of solitons by compensating the dispersion with the torque produced by 
the rotation of the dipoles induced in the medium. In both Hasegawa's theoretical model 
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(a) 




FIG. 1. (a) x-component of E\(Ef) exhibiting soliton behaviour with two maxima (twin pulses) 
when n 1R = 0.08, ftu = 0.003, r]i R = 0.02 , A = 0.1 and A = 0.3010, (b) the form of NLS 
one-soliton for comparison. 

for optical soliton in the pico-second region (temporal coherence) and in the case of 
Zhakarov/Manakov's model for modulation of the EMW in the form of solitons (spatial 
coherence), the soliton is governed by NLS equation. 

However, in our problem where the Kerr nonlinearity is not taken into account, but, by 
getting support for nonlinearity from the rotation of induced dipoles in the medium due to 
the torque with the external field the soliton is governed by the DNLS equation. The spatial 
soliton is formed when the EMW induces a waveguide in the dielectric medium via the 
nonlinearity of the torque induced, self trapped and in turn is guided in its own wave guide. 
Unlike the NLS-one soliton, here we have multihump (two) soliton (see Figs. (la) and (lb)). 
In nonlinear optics, multihump solitons corresponding to multimodes propagating without 



interference and loosing their stability are known |21]j2^]. Even though we get multihump 
(two humps) solitons, in our case it corresponds to a single mode which is again a stable 
soliton corresponding to the completely integrable DNLS equation. 



Now to understand the effect of the nonlinear term {\ip\ 2 ip)^) in the optical soliton for- 
mation during propagation of EMW we consider Eq.(^|) by dropping the dispersion term 

121- 



^ r + A(H 2 </0 ? = o. (io) 

Assuming the solution in the form of ?/>(£, r) = ?(£, T )exp[i9(£, T )}, Eq.fllOD can be written 
after rescaling r suitably as q T +3q 2 q^ = , 9 T +q 2 0^ = 0, for which the solution can 
be written in the implicit form. Infact the first of the above equations is equivalent to 
the dispersionless modified Korteweg - de Vries (mKdV) equation for which the solution 
can be written in the implicit form as q 2 = a(£ — 3rq 2 ) or equivalently in another form as 
£ = 3rq 2 + b(q 2 ), where a is an arbitrary function determined by the initial profile and b is 
the inverse function b = -. The second equation of the above equations can be solved using 
the method of characteristics and again the solution can be written in terms of another 
arbitrary function. Explicit solutions can be found for specific initial pulse profiles. For 
example, if we choose a Gaussian pulse in the form g 2 (£,0) = expi^A-) and using this, 

■so 

we find that a(£) = exp( z A-) and hence b(q 2 ) = ±£o(^i-r) 5 - This leads to the solution 
q 2 = exp[— £ 2 g ] which gives £ = 2>rq 2 ±£o(^^) 5 ; where the signs (+) and (— ) refer 
to the trailing and leading edges of the soliton respectively. From the above expression 
for £ it is clear that the local increment in velocity is inversely proportional to the square 
of the amplitude thus leading to a shift of the pulse peak to the trailing edge. As in 
the case of mKdV equation here also the self-steepening of the trailing edge leads to the 
appearance of an EM shock wave. If we repeat the calculations for a sech-profile in the form 
<Z 2 (£>0) = sech 2 (4-), the results show the formation of a shock EMW at the trailing edge of 
the pulse as in the Gaussian case. Further, if we choose a periodic wave profile in the form 
g 2 (£,0) = cos (4-), once again we obtain the same result exhibiting shocks in the negative 
slope region of the propagating EMW. This EM shock is compensated by the dispersion to 
form a localized spatially coherent optical soliton. 

In this paper we have proposed a simple and new approach to generate optical solitons in 

9 



a dielectric medium during propagation of EMW without taking into account the nonlinear 
Kerr effect as done normally. On the other hand we considered the torque produced in 
the medium due to the rotation of the induced dipoles by the interacting exterior electric 
field component of the EM field and constructed the Larmor equation for precession of the 
induced dipoles (or electric polarization) in the medium. This is then solved with Maxwell 
equations which govern the propagation of EMW in the dielectric medium using a reductive 
perturbation method. The equations were reduced to the integrable derivative nonlinear 
Schrodinger equation that admits iV-soliton solutions. Here the nonlinearity required to 
support optical soliton is provided by the torque developed between the induced electric 
dipoles in the medium and the EM field and not from the Kerr effect. Another important 
difference, when comparing the results of the usual optical soliton study, is that, unlike 
the other case here we get multihumped optical solitons of the DNLS equation. A careful 
analysis of the final results showed that the nonlinearity due to the torque actually produced 
a shock in the propagating EMW and when it is balanced by the linear dispersion, the plane 
EMW is modulated into spatially coherent optical solitons. Also, the polarization of the 
medium is excited in the form of solitons. We expect our method to bring out fascinating 
results of energy sharing optical solitons and optical logic gate operations in birefringent 
fiber medium in a straightforward way through this analytic approach. Further, this is 
expected to be useful in analysing all higher order effects in polarization of the medium. 

This work was done within the framework of the Associateship Scheme of the Abdus 
Salam International Centre for Theoretical Physics, Trieste, Italy. V.V acknowledges CSIR 
for financial support in the form of a Senior Research Fellowship. The work of M.D forms 
part of a major DST project. 
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